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ABSTRACT 

We present the latest results of our fast physical optics simulations of the ESA PLANCK HFI beams. The main 
beams of both polarized and non-polarized channels have been computed with account of broad frequency bands 
for the final design and positions of the HFI horns. Gaussian fitting parameters of the broadband beams have 
been presented. Beam polarization characteristics and horn defocusing effects have been studied. 
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1. INTRODUCTION 

With increasing the size and complexity of the space-based antennas for radioastronomy, there is a demand for 
increasing the efficiency of simulations of large multi-reflector telescope systems. A typical example is the ESA 
PLANCK Surveyor, the 3rd generation deep-space submillimeter-wave telescope being designed for measuring 
the temperature anisotropies and polarization of the Cosmic Microwave Background (CMB). The telescope 
will be equipped with two focal plane instruments, the Low-Frequency and High-Frequency Instruments, for 
detecting the radiation in the frequency range from 30 GHz to 1000 GHz in nine bands. 1 

The High-Frequency Instrument (HFI) will operate in six frequency channels 2 centered at 100, 143, 217, 353, 
545 and 857 GHz. Four channels (100 — 353 GHz) will use mono-mode quasi-Gaussian horns, 3 half of them 
feeding polarization sensitive bolometers 4 (PSB). The other two channels (545 and 857 GHz) are feeding non- 
polarized bolometers and use profiled multi-mode horns. All the HFI horns are broadband, with the bandwidth 
being about 30% of central frequency. The system is designed to meet the extreme requirements on both the 
primary mirror edge taper (—25 dB) and the angular resolution on the sky (about 5 arcminutcs at the frequencies 
of 217- 857 GHz). 

Due to the off-axis horn positions and the dual-mirror design of telescope, the beam shapes on the sky are 
far from circular (see Table 2). Their exact representation requires mapping with a significant number of points. 
This becomes important because the simulation and data reduction of the Planck mission will be huge tasks 
requiring enormous computing power. A first step towards feasibility of computations is to model the beams by 
analytical functions that can be handled more efficiently than the raw data. 5 

The aim of this paper is to present the results of the PLANCK HFI beam computations and to provide the 
update to the earlier published beam data 6 that can be used as a reference for simulations of the Planck mission 
and for the development of its data processing tools. 

2. SIMULATIONS OF THE MAIN BEAMS 

Simulations of the PLANCK HFI beams are extremely challenging because of asymmetric dual-reflector geometry 
of the telescope, large primary mirror having projected diameter D = 1.5 m (D/X m i n = 5000), a very wide field of 
view, broadband and multi-mode structure of the horn fields, and strict requirements on the accuracy. Physical 
optics (PO) is the most adequate technique for this kind of computations. Conventional software, however, 
cannot cope efficiently with problems of this size in full PO+PO mode of simulations. 
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As an alternative, we developed a dedicated fast PO simulation method which is specifically designed for large 
multi-beam multi-reflector systems with broadband channels and multi-mode structure of the horn field. 7, 8 The 
method allows us to perform rigorous PO+PO TE/TM-mode simulations of the main beams of large telescopes 
such as PLANCK in a few minutes for mono-mode mono-frequency channels (as the beams transmitted by 
telescope) and in about an hour for the broad-band polarization-averaged multi-mode channels of the highest 
frequencies. 

We compute the beam patterns of the IQUV Stokes parameters on the sky by propagating the source field 
from the apertures of corrugated horns through the telescope mode-by-mode, with integration over the frequency 
band and with account of all polarizations of the field of non-polarized channels. 8 The aperture field of horns 
has been computed by the scattering matrix approach. 9 The effective modes of the electric field at the horn 
aperture, E nm , are represented via the canonical TE-TM modes £ n j of a cylindrical waveguide as follows 

E nm (p,ip)= ^ S nm j£ nj (p,ip) (1) 

j=l,...,2M 

where S nm j — S nm j(f) is the scattering matrix provided by E. Gleeson 9 for each horn at various frequencies / 
{Snmj is used as an input in this work), n = 0, 1, N is the azimuthal index and to, j — 1,2, ...,2M are the 
radial indices accounting for both the TE (to, j = 1, M) and TM (m,j = M + 1, 2M) modes. 

The beam data are computed assuming smooth telescope mirrors of ideal elliptical shape, of perfect electrical 
conductivity of their reflective surfaces, and of ideal positioning of mirrors and horn antennas. The total beam 
power is found as a sum of powers of all modes propagated to the sky with account of all contributing polarization 
directions. Similarly, the broadband beam patterns of Stokes parameter are the sums of the relevant patterns of 
all the modes at all the frequencies of the bandwidth. In mono-mode beams, the E nm modes sum up effectively 
to a single mode which is of one unit of total power, almost Gaussian in shape and, for the polarized channels, 
of nearly perfect linear polarization on the horn aperture. 

Our latest simulations of the PLANCK HFI beams provide the updates to the beam data published in Ref. 6. 
There are five essential corrections to those data: (1) the horn design is slightly changed (mainly, due to the 
manufacturing procedure) so that the horns arc slightly elongated compared to those in Ref. 6, (2) the horn 
positions are now the final ones, (3) both broadband and monochromatic beams are computed, (4) root-mean- 
square beam fitting is performed, (5) the gains of mono-mode non-polarized beams are corrected (they should 
be increased in Ref. 6 by 3 dB while the gains of multi-mode horns remain correct for those horns at those 
positions). 

The horn positions are specified by the aperture refocus parameter Ra = Rf + Rc where Ra is the distance 
along the horn axis from the reference detector plane to the horn aperture, Rf is the similar distance to the 
point F of the geometrical focus of telescope on this axis, and Rc is the distance from the point F to the 
horn aperture. The best refocus Rao is found by minimizing the angular width of the broad-band beams 
(simultaneously, it appears to maximize the gain, even in the case of complicated multi-mode beams). At the 
best refocus, Rco = Rao — Rf specifies the focal center of the horn. This is the point inside the horn which is 
superimposed with the telescope focal point F when R A — Rao- 

The horn positions used in this paper are those of the final design. They differ in some cases from the 
ones used in Ref. 6. The final positions are specified by the values of Rc in Table 1. Generally, these Rc are 
sufficiently close to the best refocus values Rco found for the broadband horns of the final design (Fig. 1). 



Table 1. Position parameter Rc used in simulations (p or n means polarized or non-polarized horns) 



HFI beam 


100-1,4 


100-2,3 


143-p 


143-n 


217-p 


217-n 


353-p 


353-n 


545-n 


857-n 


Rc, mm 


1.2 


1.4 


1.2 


3.5 


1.1 


1.3 


1.5 


1.5 


4.0 


4.0 



Measuring polarization with bolometric detectors requires comparisons of signals of four polarization chan- 
nels obtained from the same pixel on the sky. To minimize polarization errors due to mismatch of different 
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Figure 1. HFI-545-4 horn defocusing effects: broadband beam (a) width W^b, W'so and Wf (solid, dashed and dotted 
curve, respectively), (b) gain Gb and Gf (solid and dashed curve, respectively), and (c) ellipticity sf (see Section 4) as 
a function of the horn aperture refocus parameter Ra (point shows the design value of Ra when Rc = 4.0 mm). 

beams, pairs of orthogonal polarization channels (a, b) are built into each horn by using polarization-sensitive 
bolometers. 4 The two other measurements needed to evaluate the Q and U Stokes parameters are obtained less 
than a second later, when the sister beam observes the same pixel, due to the spinning of the satellite. 

Broad frequency bands improve the power and polarization patterns of the HFI beams compared to those in 
Refs. 6: the difference of the broadband / patterns of two channels of the same beam is less than 1% of maximum 
power (Fig. 2, a). The mismatch of sister beams when superimposed on the sky is much greater (up to 5 — 8%) 
and depends on the horn location in the focal plane. 6, 8 The difference of the mono-frequency and broadband 
/ patterns of the same beam normalized to the unit total power is, typically, about 1 — 2% (Fig. 2, b). 
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Figure 2. The difference of the broadband intensity patterns of (a) two channels of orthogonal polarization of the same 
beam HFI-143-4-a/b and (b) mono- frequency and broadband beam HFI-143-4a. 

3. GAUSSIAN FITTING OF THE PLANCK HFI BEAMS 

As a development of previous work, 6 in this paper we compute the broadband HFI beams (Fig. 3). The 
broadband computations are particularly important for the multi-mode beams which have a complicated modal 
structure depending essentially on the frequency within the band. These beams are computed using up to 11 
modes in the band / = 455 — 635 GHz for the HFI-545 beams and up to 16 modes in the band / = 716 — 998 
GHz for the HFI-857 beams. All the HFI beams are computed with nine sampling frequencies in each frequency 
band that is sufficient for the broadband simulations in the limited area of the main beams. The main results 
including the Gaussian fitting parameters of the beams are presented in Table 2. 
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Figure 3. Broadband power patterns of (a) HFI-100-3a and (b) HFI-545-4 beams computed for the final design and 
positions of horns (cf respective mono-frequency patterns in Ref. 6) . 



The IQUV beam patterns are presented as functions of the (ipscVsc) coordinates on the sky {rise = 
90° — Osc) as viewed from the sky to the telescope in the spherical frame of spacecraft (SC), with the azimuthal 
and polar angles psc an d Osc, respectively, and with the polar axis being the nominal spin axis of telescope 
(the center of the focal plane corresponds to ipse = 0° and r/sc = 5°). In this representation, the Q and U 
Stokes parameters, being frame-dependent, are defined with respect to ipse as the first axis (pointing to the 
right when viewed from the sky to the telescope) and r/sc as the second axis (pointing upwards). 

The polarization angle $b of the beam field is measured from the positive direction of the (psc axis to the 
direction of major axis of polarization ellipse of the electric field E, with the positive angles counted towards the 
positive direction of the r/sc axis. So, ^e is measured counter-clockwise from the parallels of the SC frame as 
viewed from the sky in accordance with the definition of Q and U. In this case, ta,n(2^ e) = U /Q. This equation 
allows one to compute the polarization angle of the polarized component of partially polarized incoherent beams 
using the broadband values of U and Q. 

The nominal values of * B for the a-channels are ty E — 135° for the beams HFI-143-1/2, 217-5/6 and 353-3/4, 
* B = 90° for the beams HFI-143-3/4, 217-7/8 and 353-5/6, and * B = 112.5°, 135°, 90° and 67.5° for the beams 
HFI-100-1, 2, 3 and 4, respectively. The complementary b-channels have the PSB polarization directions in the 
horns precisely orthogonal to those of a-channcls. On the sky non-orthogonality of polarization is less than 
0.03° for the beam-average angles * E and, typically, less than 0.10° for the polarization angles on the beam 
axes. Generally, the beam-average angles ^e coincide with the required nominal values better than by 0.01°, 
though some on- axis values may differ by about 0.1°. 

For the polarized beams, along with the distortion of polarization, some depolarization appears when the 
field propagates from the ideal PSB through the horn and the telescope. Similarly, for the non-polarized beams, 
small polarization arises. These effects are represented by the deviations of the QUV Stokes parameters from 
the ideal values across the beam patterns. The bounds on these non-idealities are found as 5V = max(\V\/ I max ) 
and SL = max((5Q, SU) where SQ = max(|Q — Qq\/ I max ), SU — max(|£/ — Uo\/I max ) and I m ax = max(J), with 
Qo and Uq being the ideal values of Q and U, respectively (the values when the degree of polarization is either 
one or zero and the polarization angle ^e is equal to its nominal value across the whole beam). Notice, that 
the beam-average non-idealities SVb and 5 Lb are much smaller than the maximum values SV and 5L. 

In this work, we fit the intensity Stokes parameter patterns I(psc,Vsc) by the elliptical Gaussian beams 
F(<Psc,Vsc)- We normalize all the Stokes parameters so that max(7) = 1. In the original PO simulated beams, 
the power per one steradian, P [Watt/sr], is computed with account of all contributing modes and polarizations, 
in the units of power Pq [Watt / sr] radiated per one steradian by an isotropic source with the power of one mode 
of single polarization. Then, the original patterns P(psc, Vsc), when presented in dBi, are computed as follows 



P(¥sc,Vsc)[dBi] = 101ogi {/ (ipse ,Vsc)} + ^WdBi] (2) 

where P ma x [dBi] is the power gain of the original beam in the units of power P of an isotropic source of one 
mode of single polarization. The power of one mode constituting the unit Po is defined as being radiated (or 
received) by the bolometer of the given cross-section sb accounted in the scattering matrix S nm j. 

With the incident radiation characterized by the Stokes parameters {I s ky, Qsky, U s k v , V s k y } (in absolute units, 
so that I s fc y [Watt/(sr • Hz)] is the sky brightness), the total power S a [Watt] absorbed by the bolometer in the 
frequency band A/ (the channel response) is evaluated as 



S a = 0.5P max j df J dVL{II sky + QQ s k v + UU sky + VV sky } 



where P max [rel.un.] is related to the beam gain P max [dBi] in Eq. (2) as P max [dBi] = 101ogi {P„ 
and d£l is the element of solid angle (the gains in Table 2 are given per unit frequency range) . 

For the beam fitting, we use the elliptical Gaussian function 

F&scVsc) = Aexp[-(p 2 /a 2 + t 2 /b 2 )}, 



(3) 

.[rel.un.]} 



(4) 



where p = xcos(r) + ys m(r), t = -xsin(r) + ycos(r), x = ipse - <^sco, V = Vsc - Vsco, a = 0.5W max /q, 
b = 0.5W m i n /q, q = y/ln(2), (psco an d rjsco are the angular coordinates of the center point (the point of 
maximum power) of the fitting elliptical Gaussian beam in the SC frame, W max and W m i n are the full beam 
widths at half magnitude measured along the major and minor axes of the fitting beam ellipse, respectively, and 
r is the angle from the ipse ax is to the major axis of the beam ellipse measured counter-clockwise as viewed 
from the sky to the telescope (Fig. 4, a). 

Fit error: HFI-100-3-a. f= 84 ..116 GHz 
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Figure 4. (a) Parameters of the elliptical fitting beam and (b) fitting error of the broadband beam HFI-1003-a (cf Ref. 6). 

According to this definition, for each HFI beam, there are six fitting parameters (W max , W m i n , ipsco, Vsco, 
t, and A) which can be found by minimizing the relevant aim function that quantifies the deviation of the 
elliptical Gaussian fit from the given PO simulated beam. One can propose a variety of the aim functions to 
evaluate the difference between the fit and the actual beam. Because the basic quantity of interest is the power 
contribution to the PSB readout, it is the deviation of the beam power profile from the Gaussian fit that has to 
be minimized. Still, even in this special case, the aim function can be defined in various ways, see Ref. 6. 

Unlike the options listed in Ref. 6, in this paper we use the root-mean-square (RMS) misfit of the beam 
pattern as the aim function. The latter is defined as 

AF= ( f(F- Ifdfl ) / I [ IdSl, (5) 



though the fitting error in Table 2 is estimated as the maximum of SF = \F — I\. In addition, we require the 
conservation of the total power of the beam Pb = J FdD, = J Id£l so that the number of independent fitting 
parameters reduces to five. 

4. RESULTS AND COMMENTS 

The RMS Gaussian fitting parameters of the broadband HFI beams are shown in Table 2. Table 2 shows also 
some other parameters of the HFI beams such as the full width at the 3dB level of the original PO simulated beam 
(Wb = Wz<ib), the mean full width at half magnitude of the Gaussian fitting beam Wf — ^Wmax^min) 1 ^ 2 , 
the original beam gain Gs[dBi] = P TOOa; [dBi], the fitting beam gain GV[dBi] = P max [dBi] + 101og 10 ( J 4), the 
fitting beam ellipticity ep = W max /W m i n , the maximum fit error SF, and the bounds SL and SV on the peak 
deviations of the QUV Stokes parameters of the original beams from the ideal values. 

Overview of the data in Table 2 shows that the peak errors of the elliptical Gaussian fit with respect to the 
original beams is, typically, about 2% of maximum power for the mono-mode quasi- Gaussian beams HFI-100, 
143, 217 and 353 (Fig. 4, b), while being up to 40% for the essentially flat-top multi-mode beams HFI-545 
and 857. The error of 2% is quite significant, being comparable to the difference between mono-frequency and 
broad-band beams of polarized channels and usually greater than the difference between the beams of orthogonal 
polarizations of the same horns (see Section 2). This means that the elliptical Gaussian fit should be considered 
as a first term of a more advanced fitting to be developed later. 

On the other hand, the fitting error remains smaller than the typical difference between the beams of two 
different complementary horns (the latter is about 5%). It means that the elliptical Gaussian fits are quite 
capable of accounting for the beam mismatch effects essential for polarization measurements. 

Another essential observation is that the beam width of the Gaussian fit, Wf, is noticeably smaller than the 
typical width of the original beam, Wb , evaluated as a diameter of the circle of the same area as bounded by the 
3dB isolevel (Wb — Wz<ib), though it is close to the diameter of the area bounded by the isolevel encompassing 
50% of the total power of the beam (Wb — VF50), see Fig. 1, a. The reason is that the real beams have 
specific shape at the top, tending to be flat-top, though more complicated. This happens even to quasi-Gaussian 
beams (more so to the broadband ones, see Fig. 2, b) both due to the field propagation via the horns and 
via the telescope, though multi-mode beams are significantly more flat-top because of their complicated modal 
composition. As a result, the original beams at half magnitude are always wider compared to the Gaussian fits 
at their half magnitude, respectively. More complicated beam fitting, when using additional terms, would be 
able to represent the width of the original beams with a better accuracy. 

5. CONCLUSION 

We present the results of our simulations of the broadband HFI beams of the ESA PLANCK Surveyor. The beams 
are computed by multi-mode physical optics propagation of the source field from the apertures of corrugated 
horns simulated by the scattering matrix approach. The horns of the final design at the actual positions in the 
focal plane of the telescope have been used. 

Broad frequency bands of the HFI channels are shown to improve both the beam power and polarization 
patterns. Perfect alignment of polarization along the required directions is confirmed by small and symmetric 
deviations of the respective Q or U patterns from the ideal values. Peak difference of power patterns of orthogonal 
polarization channels, due to the usage of orthogonal pairs of polarization sensitive bolometers in the individual 
horns, is less than 1%, while it is 5% to 8% for the different beams of complementary pairs of horns designed 
for the polarization measurements. 

Compact representations of the IQUV patterns by a few parameters of the elliptical Gaussian fitting function 
are found. This approximation is generally considered as acceptable from the scientific viewpoint, although we 
show that induced errors are far from negligible. 

New simulation technique proves indispensable for multi-reflector systems with a very wide field of view, 
large number of beams and complicated structure of propagating waves, when efficient and rigorous simulations 
of the main beams are needed. 
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Table 2. Gaussian fitting parameters and some representative characteristics of the broadband PLANCK HFI beams 
(horns are located at the final design positions specified with respect to the reference detector plane, beams are computed 
with nine sampling frequencies, Wb = W 3d B, Wf = {W m axW m i n ) 1/2 ) 
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SL 


8V 
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arcmin 


arcmin 


dBi 
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% 


% 
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100-la 


10.4560 


8.9481 


1.1842 


3.5120 


13.5124 


9.90 


9.67 


60.33 


60.40 


1.169 


1.9 


0.4 


2.1 


100- lb 
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8.9528 


1.1842 
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13.4052 


9.88 


9.67 


60.33 
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1.8 


0.4 


2.0 


100-2a 


10.5047 


8.9918 


0.3829 


3.2410 


4.4241 


10.05 


9.72 


60.61 


60.70 


1.168 


2.1 


0.4 


1.3 


100-2b 
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3.6527 
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100-3a 


10.5147 
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1.3 


100-3b 
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10.05 
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1.3 


100-4a 


10.4579 
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-1.1842 


3.5120 
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2.1 


100-4b 


10.4525 
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8L 
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8V 
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143-2a 

143-2b 

143-3a 

143-3b 

143-4a 

143-4b 

143-5 

143-6 

143-7 
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7.4323 
7.4986 
7.1939 
7.2439 
7.2249 
7.1967 
7.5149 
7.5083 
7.5902 
7.3981 
7.3836 
7.6067 



6.8745 
6.8574 
6.9619 
6.9324 
6.9377 
6.9460 
6.8715 
6.8790 
7.1702 
7.1650 
7.1576 
7.1833 



1.3718 
1.3717 
0.5637 
0.5637 
-0.5637 
-0.5637 
-1.4418 
-1.4419 
1.1443 
0.2997 
-0.2997 
-1.1443 



6.1966 
6.1966 
6.2243 
6.2243 
6.1992 
6.1993 
6.2210 
6.2211 
6.7334 
6.7608 
6.7358 
6.7585 



49.6927 
49.0453 
61.3503 
56.6716 
-61.7208 
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-51.4655 
-47.3673 
62.8153 
79.0498 
-78.5876 
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7.30 
7.35 
7.29 
7.32 
7.32 
7.29 
7.33 
7.33 
7.50 
7.41 
7.41 
7.51 



7.15 
7.17 
7.08 
7.09 
7.08 
7.07 
7.19 
7.19 
7.38 
7.28 
7.27 
7.39 



63.44 
63.41 
63.54 
63.52 
63.53 
63.54 
63.40 
63.40 
66.23 
66.34 
66.35 
66.21 



63.51 
63.48 
63.62 
63.60 
63.61 
63.62 
63.47 
63.47 
66.29 
66.41 
66.42 
66.27 



1.081 
1.094 
1.033 
1.045 
1.041 
1.036 
1.094 
1.091 
1.059 
1.033 
1.032 
1.059 



2.1 
2.0 
2.2 
2.2 
2.2 
2.2 
2.0 
2.0 
2.1 
2.0 
2.0 
2.1 



0.9 
1.0 
0.6 
0.5 
0.7 
0.5 
1.2 
1.0 
0.5 
0.5 
0.5 
0.5 



4.0 
4.2 
3.7 
3.8 
3.7 
3.8 
4.2 
4.3 
4.5 
4.2 
4.2 
4.5 



217-1 

217-2 

217-3 

217-4 

217-5a 

217-5b 

217-6a 

217-6b 

217-7a 

217-7b 

217-8a 

217-8b 



5.0748 
5.0071 
5.0086 
5.0725 
5.0741 
5.0718 
4.9771 
4.9756 
4.9781 
4.9755 
5.0733 
5.0716 



4.4924 
4.5114 
4.5070 
4.4947 
4.5330 
4.5339 
4.5835 
4.5841 
4.5773 
4.5798 
4.5354 
4.5377 



0.9960 
0.3149 
-0.3148 
-0.9961 
1.2256 
1.2257 
0.5435 
0.5436 
-0.5435 
-0.5435 
-1.2258 
-1.2258 



4.0142 
4.0403 
4.0153 
4.0391 
4.5149 
4.5149 
4.5416 
4.5416 
4.5166 
4.5166 
4.5399 
4.5399 



12.1172 
4.3963 
-4.3277 
-12.3129 
18.5837 
18.3535 
10.9266 
10.5525 
-10.4282 
-10.5638 
-18.7074 
-18.8780 



4.88 
4.86 
4.85 
4.88 
4.90 
4.90 
4.87 
4.87 
4.87 
4.88 
4.90 
4.90 



4.77 
4.75 
4.75 
4.77 
4.80 
4.80 
4.78 
4.78 
4.77 
4.77 
4.80 
4.80 



70.03 
70.08 
70.08 
70.03 
66.96 
66.97 
67.01 
67.01 
67.02 
67.02 
66.97 
66.96 



70.10 
70.15 
70.15 
70.10 
67.02 
67.03 
67.08 
67.08 
67.09 
67.09 
67.03 
67.03 



1.130 
1.110 
1.111 
1.129 
1.119 
1.119 
1.086 
1.085 
1.088 
1.086 
1.119 
1.118 



2.0 
2.1 
2.1 
2.0 
2.2 
2.2 
2.2 
2.2 
2.2 
2.2 
2.2 
2.2 



0.2 
0.2 
0.2 
0.2 
0.4 
0.4 
0.3 
0.3 
0.3 
0.2 
0.4 
0.4 



1.5 
0.5 
0.5 
1.5 
2.0 
2.0 
1.2 
1.2 
1.2 
1.2 
2.0 
2.1 



353-1 

353-2 

353-3a 

353-3b 

353-4a 

353-4b 

353-5a 

353-5b 

353-6a 

353-6b 

353-7 

353-8 



5.1917 
4.9061 
4.7513 
4.7513 
4.6682 
4.6681 
4.6832 
4.6822 
4.7825 
4.7812 
4.9395 
5.1953 



4.1539 
4.2351 
4.3493 
4.3500 
4.4435 
4.4437 
4.4224 
4.4220 
4.3264 
4.3252 
4.2137 
4.1554 



2.0626 
1.4195 
0.8248 
0.8248 
0.2107 
0.2108 
-0.3687 
-0.3687 
-0.9650 
-0.9650 
-1.5240 
-2.0628 



5.0155 
5.0436 
5.0204 
5.0204 
5.0462 
5.0462 
5.0211 
5.0211 
5.0451 
5.0451 
5.0181 
5.0406 



21.7076 
21.9361 
19.6160 
19.5862 
8.4538 
8.3505 
-12.6675 
-12.6907 
-21.0534 
-21.0458 
-21.4681 
-22.1009 



4.66 
4.66 
4.66 
4.66 
4.64 
4.64 
4.64 
4.64 
4.65 
4.64 
4.64 
4.69 



4.64 
4.56 
4.55 
4.55 
4.55 
4.55 
4.55 
4.55 
4.55 
4.55 
4.56 
4.65 



70.58 
70.74 
67.68 
67.68 
67.67 
67.67 
67.68 
67.68 
67.68 
67.68 
70.73 
70.58 



70.60 
70.79 
67.74 
67.74 
67.74 
67.74 
67.75 
67.75 
67.75 
67.75 
70.78 
70.61 



1.250 
1.158 
1.092 
1.092 
1.051 
1.051 
1.059 
1.059 
1.105 
1.105 
1.172 
1.250 



2.0 
1.6 
1.6 
1.6 
1.6 
1.6 
1.6 
1.6 
1.6 
1.6 
1.7 
2.0 



0.1 
0.1 
0.7 
0.6 
0.5 
0.4 
0.3 
0.4 
0.5 
0.5 
0.1 
0.1 



2.4 
1.8 
1.3 
1.3 
1.2 
1.2 
1.2 
1.1 
1.6 
1.5 
1.9 
2.4 



545-1 
545-2 
545-3 
545-4 



4.4859 
3.8050 
3.8636 
4.5011 



3.7289 
3.3962 
3.4156 
3.7395 



2.0686 
1.4230 
-1.5279 
-2.0688 



5.5336 
5.5632 
5.5376 
5.5585 



35.5850 
41.3291 
-40.2154 
-35.9540 



4.77 
4.40 
4.40 
4.82 



4.09 
3.59 
3.63 
4.10 



75.74 
76.84 
76.75 
75.71 



76.33 
77.47 
77.38 
76.30 



1.203 
1.120 
1.131 
1.204 



25.8 
40.3 
38.6 
25.6 



1.9 
3.1 
2.9 
1.9 



7.0 
7.1 
7.3 
7.0 



857-1 
857-2 
857-3 
857-4 



3.9595 
3.8445 
3.8642 
3.9924 



3.6561 
3.7628 
3.7386 
3.6410 



0.8440 
0.2288 
-0.3519 
-0.9669 



5.5381 
5.5641 
5.5390 
5.5628 



37.0723 
37.9611 
-37.0791 
-37.4161 



4.93 
4.95 
4.91 
4.93 



3.80 
3.80 
3.80 
3.81 



79.77 
79.73 
79.75 
79.73 



80.46 
80.47 
80.48 
80.43 



1.083 
1.022 
1.034 
1.097 



36.6 
37.3 
37.3 
36.7 



10.0 2.7 

10.3 2.2 

10.3 2.1 

9.8 3.2 



